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Exclusive photoproduction of dileptons, 7A'' — > l^l" N, is and will be measured in ultraperipheral 
collisions at hadron colliders, such as the Tevatron, RHIC and the LHC . We demonstrate that 
the timelike deeply virtual Compton scattering (TCS) mechanism 75 £'^£~q where the lepton 
pair comes from the subprocess 7^ j*q dominates in some accessible kinematical regions, thus 
opening a new way to study generalized parton distributions (GPD) in the nucleon. High energy 
" kinematics enables to probe parton distributions at small skewedness. This subprocess interferes at 

. the amplitude level with the pure QED subprocess 77* — > £^£~ where the virtual photon is radiated 

^ ' from the nucleon. 

^ ' PACS numbers: 13.60.Fz , 13.90.+i 
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! I. INTRODUCTION. 



f~i ■ Much theoretical and experimental progress has recently been witnessed in the study of deeply virtual Compton 
CLh; scattering (DVCS), i.e., 7*p "fp, an exclusive reaction where generalized parton distributions (GPDs) factorize 
from perturbatively calculable coefficient functions, when the virtuality of the incoming photon is high enough 1]. It 
is now recognized that the measurement of GPDs should contribute in a decisive way to our understanding of how 
, quarks and gluons build hadrons [3] . In particular the transverse location of quarks and gluons become experimentally 
■ measurable via the transverse momentum dependence of the GPDs [3| . 
The "inverse" process, 

-f{q)N{p) r{q')N{p') -> l-{k)l+{k')N{p') 

I at small t = {p' — p)'^ and large timelike virtuality (fc + fc')^ = q'^ ~ Q'^ of the final state dilepton, timelike Compton 
i scattering (TCS) shares many features with DVCS. The Bjorken variable in that case is r = Q'^/s with s = (p+q)^- 
One also defines A = p' — p {t = A'^) and the skewness variables ^ , 77 as 
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5h ' where the approximations hold in the kinematical limit we are working, i.e. in the extended Bjorken regime where 



masses and —t are small with respect of Q'^ (s is always larger than Q'^ ). x, ^, and r] represent plus-momentum 
fractions (Light-cone coordinates are defined as = ^ , both proton momenta p and p' moving fast to the right, 
i.e., having large plus-components). 

(p + p')+ {P + P')+ {P + P')+' 

To leading-twist accuracy one has f = —77 = —t/{2 — r). 

The possibility to use high energy hadron colliders as powerful sources of quasi real photons in ultraperipheral 
collisions has recently been emphasized Q ■ This should allow the study of many aspects of photon proton and photon 
photon collisions at high energies, already at the Tevatron and at RHIC but in particular at the LHC ^ even if 
the nominal luminosity is not achieved during its first years of operation. The high luminosity and energies of these 
photon beams opens a new kinematical domain for the study of TCS , and thus to the hope of determining GPDs in 
the small skewedness region, which is complementary to the determination of the large ^ quark GPDs at lower 
energy electron accelerators such as JLab. Moreover, the crossing from a spacclike to a timelike probe is an important 
test of the understanding of QCD corrections, as shown by the history of the understanding of the Drell-Yan reaction 
in terms of QCD. 



Figure 1: Real photon-proton scattering into a lepton pair and a proton. 




Figure 2: Kinematical variables and coordinate axes in the 7p and i cm. frames. 



The physical process where to observe TCS is photoproduction of a heavy lepton pair, — )■ fi^fi^ N or 'jN — > 
e+e" N, shown in Fig. [TJ As in the case of DVCS, a Bethe-Heitler (BH) mechanism - sometimes called 77 process 
since the lepton pair is produced through the '-f{q)'-f{A) subprocess - contributes at the amplitude level. This 

amplitude is completely calculable in QED provided one knows the Nucleon form factors at small t. This process 
has a very peculiar angular dependence and overdominates the TCS process if one blindly integrates over the final 
phase space. One may however choose kinematics where the amplitudes of the two processes are of the same order 
of magnitude, and either subtract the well-known Bethe-Heitler process or use specific observables sensitive to the 
interference of the two amplitudes. 

The kinematics of the ^{q)N{p) £^ {k)£^ {k')N {p') process is shown in Fig. [2l In the £^£^ center of mass system, 
one introduces the polar and azimuthal angles and ip of fc, with reference to a coordinate system with 3-axis along 
—p' and 1- and 2-axes such that p lies in the 1-3 plane and has a positive 1-component. 

In this paper, we shall examine in detail the feasibility of TCS experiments in ultraperipheral collisions at the LHC. 
Most conclusions should also apply to the Tevatron and RHIC experimental conditions, once the effective photon 
fluxes and energies are scaled down to their specific values. We shall work in the leading twist approximation. As in 
the DVCS case, a gauge invariant treatment necessitates the inclusion of twist 3 effects J,], but we will not address 
this problem in this paper. We will also stay at the leading order in as, and thus neglect contributions proportionnal 
to gluon GPDs. The motivation of this study has been presented at a recent LHC workshop Q. 
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Figure 3: The Feynman diagrams for the Bethe-Heitler amplitude. 

II. THE VARIOUS CONTRIBUTIONS 
A. The Bethe-Heitler contribution 

The Bethe-Heitler amplitude is calculated from the two Feynman diagrams in Fig. [3] where the photon-nucleon 
vertex is parameterized by the usual Dirac and Pauli form factors Fi{t) and F2{t), normalizing ^2(0) to be the 
anomalous magnetic moment of the target. Neglecting masses and t compared to terms going with s or Q'^, the 
Bcthc Heitler contribution to the unpolarized jp cross section is {M is the proton mass) 
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provided we stay away from the kinematical region where the product of lepton propagators goes to zero at very small 
9. The interesting physics program thus imposes a cut on 9 to stay away from the region where the Bethe Heitler 
cross section becomes extremely large. 



B. The Compton ampHtude 

In the region where the final photon virtuality is large, the amplitude is given by the convolution of hard scattering 
coefficients, calculable in perturbation theory, and generalized parton distributions, which describe the nonperturba- 
tive physics of the process. To leading order in as one then has the dominance of the quark handbag diagrams of 
Fig. m The analysis of these handbag diagrams show the simple relations 
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between the helicity amplitudes for TCS and DVCS at equal values of 77 and t. For instance, 
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where the Compton form factors TLi,£i,TLi,£i are defined as 
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Figure 5: The NLO(MS') GRVGJR 2008 parametrization of u{x) + u{x) for different factorization scales fi% — 4 (dotted) , 5 
(dasiied) , 6 (dash-dotted) , 10 (solid) GeV^ 



and H'^{x, ij, t), E^{x, rj, t), H'^{x, rj, t), E'^{x, rj, t) are usual GPDs for a quark of flavour q and electric charge e e^. 

At this order the two processes carry the same information on the generalized quark distributions. This will not be 
true when higher order contributions are included but we will not consider these corrections here. 

The crucial ingredient to estimate the TCS amplitude at large energies is a realistic model of GPDs at small 
skewedness. 



C. Modelizing GPDs 



We anticipate that singlet quark GPDs give the dominant contributions to the TCS amplitude in that domain. 
Since gluon GPDs only enter the TCS amplitude at the 0{as) level, we feel justified in a first step to neglect their 
contributions and leave its study for a future work. The choice of factorization scale is a major issue since GPD 
evolution is particularly active in the small x domain. In this first study of the feasibility of the extraction of the TCS 
signal, we feel justified to simplify our calculations by using a factorization ansatz for the t dependence of GPD's: 



H^X,7J,t) - h^x,Tj)-F^it) 

H^{x,V,t) = h''ix,rj)Ftit) 
H\x,r^,t) = h\x,ri)FD{t) 
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Figure 6: h\{x,ri) — h'"{x,rj) — h^{—x,ri) for i) — IQ ^ (a) and for r] — IQ ^ (b) for different factorization scales — 4 
(dotted) , 5 (dashed) , 6 (solid) GeV^ 




Figure 7: (a) The BH cross section integrated over 9 G [7r/4, 37r/4], Lp € [0,2tv] , Q'^ G [4.5, 5.5] GeV^ \t\ G [0.05, 0.25] GeV^ 
as a function of 7p cm. energy squared s. (b) The BH cross section integrated over (p G [0,27r] , \t\ G [0.05, 0.25] GeV^ and 
various ranges of e : [tt/S, 27r/3] (dotted), [7r/4, 37r/4] (dashed) and [tt/G, 5tt/6] (solid), as a function of Q' for s = 10^ GeV^ 



and a double distribution ansatz for h'' without any D-term: 



h'^{x,ri) = / dx' 
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For the unpolarized distributions q{x) and q{x) we take NLO(MS) GRVGJR 2008 parametrization [9]. Their strong 
dependence of the factorization scale choice for small x is shown on FigO This results in the strong dependence of 
/i' for small values of rj as shown on Fig [HI 



III. CROSS SECTION ESTIMATES 



Let us now estimate the different contributions to the lepton pair cross section for ultraperipheral collisions at the 
LHC. Since the cross sections decrease rapidly with Q'^, we are interested in the kinematics of moderate Q'^, say a 
few GeV^, and large energy, thus very small values of rj. Note however that for a given proton energy the photon flux 
is higher at smaller photon energy. 
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A. The Bethe Heitler cross section 



The full Bethe Heitler cross section integrated over 9 e [7r/4, 37r/4], ip e [0,27r] , Q'^ e [4.5, 5.5] GeV^, \t\ e 
[0.05,0.25] GeV^, as a function of 7p energy squared s is shown on Fig. [7^. We see that in the limit of large s 
it is constant and equals 28.4 pb. On Fig. (TJj, the Bethe Heitler contribution is shown as a function of Q'^ when 
it is integrated over (p in the range [0, 27r], —t in the range [0.05, 0.25]GeV^ and for 9 integrated in various ranges 
[tt/S, 2tt/3], [tt/A, 3tt/A] and [tt/6, 5tt/6]. As anticipated, the cross section grows much when small 9 angles are allowed. 
In the following we will use the limits [7r/4, 37r/4] where the cross section is sufficiently big and does not dominate too 
much over the Compton process. 



B. The TCS cross section 



Since the u-quark contribution of Hi turns out to dominate the TCS amplitude, we show in Fig. [5] the real and 
imaginary part of divided by ^F'^{t) for various factorization scales and various ranges of ^. We observe that for 
small values of ^, factorization scale dependence is quite strong. The same is seen on Fig. [9] where the full Compton 
cross section arcs is plotted as a function of the photon-proton energy squared s. For very high energies arcs 
calculated with = 6GeV^ is much bigger then with — 4GeV^. Also predictions obtained using LO and NLO 
GRVGJR2008 PDFs differ significantly. 



C. The interference cross section 



Since the amplitudes for the Compton and Bethe-Heitler processes transform with opposite signs under reversal of 
the lepton charge, the interference term between TCS and BH is odd under exchange of the £^ and £~ momenta. It is 
thus possible to project out the interference term through a clever use of the angular distribution of the lepton pair. 
The interference part of the cross-section for 7p — > i^i" p with unpolarized protons and photons is given at leading 
order by 
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where —to — Arf'M'^/{\ — rf) . With the integration limits symmetric about 9 — 7r/2 the interference term is odd 
under ip n + ip due to charge conjugation, whereas the TCS and BH cross sections are even. One may thus extract 



the Compton amplitude through a study of J cos < 



In Fig. [TO]we show the interference contribution to the cross section in comparison to the Bethe Heitler and Compton 
processes, for various values of photon proton energy squared s — 10^ GeV^, 10^ GeV^, 10'^ GeV^. We observe that for 
larger energies the Compton process dominates, whereas for s = 10^ GeV" all contributions are comparable. 



D. Rate estimates 



As described in 10] the cross section for photoproduction in hadron collisions is given by: 
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(9) 



where a~fp{k) is the cross section for the 7p — *■ pl~^l process and k is the photon energy, '^^^jf' is an equivalent photon 
flux (the number of photons with energy fc), and is given by [ll|: 
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Figure 8: ImTC^ (left) and RcH" (right) divided by ^F" for various factorization scales = 4 (dotted), 5 (dashed), 6 (solid) 
GeV^ and various ranges of C : [1 ■ 10"\ 2 • 10"^], [1 • 10"^, 2 • 10"-'],[1 • 10-^ 2 • 10"^]. 



where: A = 1 + ^'"^qP^^ , Qmin ~ — — minimal squared fourmomentum transfer for the reaction, and Spp is 

the proton-proton energy squared {^/s^ = 14TeV). The relationship between energy squared s and k is given by: 

s w 2^Sppk (11) 

The pure Bethe - Heitler contribution to app, integrated over 6 = [7r/4, 37r/4], ^ = [0, 27r], t = 
[-0.05 GeV^ -0.25 GeV% Q'^ = [4.5 GeV^ 5.5 GeV^], and photon energies k = [20, 900] GeV gives: 

af/ = 2.9pb. (12) 

The Compton contribution (calculated with NLO GRVGJR2008 PDFs, and n% = 5GeV^) gives: 

aj_p^ = 1.9pb. (13) 
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Figure 9: arcs as a function of -yp cm. energy squared s, for GRVGJR2008 LO (a) and NLO (b) parametrizations, for different 
factorization scales /i^ = 4 (dotted), 5 (dashed), 6 (solid) GeV^. 



We have choosen the range of photon energies in accordance with expected capabihties to tag photonenergies at the 
LHC. This amounts to a large rate of order of 10^ events/year at the LHC with its nominal luminosity (10'^'* cm~^s~^). 
The rate remains sizeable for the lower luminosity which can be achieved in the first months of run. 

IV. CONCLUSION 

The operation of LHC as a heavy ion collider will enable us to study TCS on nuclei. Such scattering may occur 
in a coherent way and its amplitude involves nuclear GPDs This very interesting subject definitely needs more 
work. Incoherent TCS which occurs on quasi free neutrons and protons allows to study the GPDs of bound nucleons. 
Cross sections arc then roughly multiplied by A^, where A is the atomic number. One may thus expect sizeable rates 
in ion-ion collisions at the LHC. 

In conclusion, let us stress that timelike Compton scattering in ultraperipheral collisions at hadron colliders opens 
a new way to measure generalized parton distributions. Our work has to be supplemented by studies of higher order 
contributions which will involve the gluon GPDs; they will hopefully lead to a weaker factorization scale dependence 
of the amplitudes. 
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